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Abstract 

The relativistic Coulomb resummation factor suggested by I.L. Solovtsov 
is used to reanalyze the mass limits obtained for magnetic monopoles 
which might have been produced at the Fermilab Tevatron. The limits 
given by the Oklahoma experiment (Fermilab E882) are pushed close 
to the unitary bounds, so that the lower limits on monopole masses are 
increased from around 250 GeV to about 400 GeV. 



1 Threshold resummation 

In describing a charged particle-antiparticle system near threshold, it is well 
known from QED that the so-called Coulomb resummation factor plays an 
important role. This resummation, performed on the basis of the nonrel- 
ativistic Schrodinger equation with the Coulomb potential V{r) — —a/r^ 
leads to the Sommerfeld-Sakharov S^-factor [1]. In the threshold region one 
cannot truncate the perturbative series and the S'-factor should be taken into 
account in its entirety. The S^-factor appears in the parametrization of the 
imaginary part of the quark current correlator, which can be approximated 
by the Bethe-Salpeter amplitude of the two charged particles, Xbs(^ — 0) 
[2]. The nonrelativistic replacement of this amplitude by the wave function, 
which obeys the Schrodinger equation with the Coulomb potential, leads 
to the appearance of the resummation factor in the parametrization of the 
normalized e^e~ to hadrons cross section ratio R{s). 

For a systematic relativistic analysis of quark-antiquark systems, it is 
essential from the very beginning to have a relativistic generalization of the 
^-factor. A new form for this relativistic factor in the case of QCD was 
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proposed by Solovtsov in 2001 [3]. 

Six) = , ^l^Vr M ' ^(^) = -^' « 

1 - exp [-X(x)] smhx 

where x is the rapidity which is related to s by 2m cosh x = <^ ^ 4^5/3 
in QCD. The function X{x) can be expressed in terms of v = — Am^/s: 
^{x) — Tvay/l — jv. The relativistic resummation factor (1) reproduces 
both the expected nonrelativistic and ultrarelativistic hmits and corresponds 
to a QCD-hke Coulomb potential. Here we consider the vector channel for 
which a threshold resummation S-f actor for the s-wave states is used. For the 
axial-vector channel the P-factor is required. The corresponding relativistic 
factor has been found recently [4]. 

To incorporate the quark mass effects one usually uses the approximate 
expression above the quark-antiquark threshold [5] 

'R{s)=T{v) \\^g{v)Tis)\ , (2) 

where 
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The function g{y) is taken in the Schwinger approximation [1]. 

One cannot directly use the perturbative expression for t{^s) in Eq. (2), 
which contains unphysical singularities, to calculate, for example, the Adler 
D-function. Instead, one can use the analytic perturbation theory (APT) 
representation for t{s). The explicit three- loop form for vp^^^is) can be 
found in Ref. [6]. Besides this replacement, one has to modify the expression 
(2) in such a way as to take into account summation of an arbitrary number 
of threshold singularities. Including the threshold resummation factor (1) 
leads to the following modification of the expression (2) for a particular 
quark flavor / [6, 7] 

7^/(5) = \R^j{s) + Rrj{s)\ 9(5 - 4m^), (4) 
i?o(5) =T(i;)5(x), RM-T[v) 



The usage of the resummation factor (1) reflects the assumption that the 
coupling is taken in the V renormalization scheme. To avoid double count- 
ing, the function R\ contains the subtraction of X{x). The potential term 
corresponding to the i?o function gives the principal contribution to 7^(5), 
as shown in Fig. 1, the correction R\ amounting to less than twenty percent 
over the whole energy interval. For a recent account of some of the successes 
of APT including Coulomb resummation see Ref. [8]. 
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Fig. 1: Behavior of Ro{s) with relativistic and nonrelativistic S -factors. 
2 Dirac Magnetic Monopoles 

The relativistic interaction between an electric and a magnetic current is [9] 

W{j, *j) = J {dx){dx'){dx"rf{x)e^,^rdT {x - x') D+ {x' - x") f {x") . 

(5) 

Here the electric and magnetic currents are 

j/x = e^p-fj^^i/j and = gxi/^X^ (6) 

for example, for spin-1/2 particles. The photon propagator is denoted by 
— x^) and fi_i{x) is the Dirac string function which satisfies the differ- 
ential equation 

d^rix) = 4irS{x). (7) 
A formal solution of this equation is given by 

/^(x) =47Tn^{n-dy^6{x), (8) 

where is an arbitrary constant vector. 

Dirac showed in 1931 [10] that quantum mechanics was consistent with 
the existence of magnetic monopoles provided the quantization condition 
holds, 

eg = m'fic, (9) 



3 



where m' is an integer or an integer plus 1/2, which explains the quantization 
of electric charge. This was generalized by Schwinger to dyons, particles 
carrying both electric charge and magnetic charge Qa [11]: 

ei^2 - 62^1 = -m'Tic. (10) 

(Schwinger sometimes argued that m' was an integer, or perhaps an even 
integer.) For details on the derivation of these quantization conditions, see, 
for example, Ref. [9]. In the following we write m! = n/2. 

3 OU Monopole Experiment: Fermilab E882 

We now refer to the experiment conducted at the University of Oklahoma 
from 1997-2004 [12], searching for low-mass monopoles which might have 
been produced at the Tevatron and captured in the old CDF and DO detec- 
tors. 

The best prior experimental limit on the direct accelerator production of 
magnetic monopoles is that of Bertani et al. in 1990 [13] 

cr < 2 X lO'^'^cm^ for a monopole mass M < 850 GeV. (11) 

(Such limits are complementary to searches for cosmic "intermediate mass" 
magnetic monopoles, with masses between 10^ and 10^^ GeV, such as have 
been recently reported in Ref. [14].) We are able to set much better limits 
than Bertani et al. because the integrated luminosity is 10^ times that of the 
previous 1990 experiment: 

y*£ = 172±8pb-^ (DO). (12) 

The fundamental mechanism is supposed to be a Drell-Yan process, 

p + p-^q + q + X-^M + M + X, (13) 
where the cross section is given by 

Here M is the invariant mass of the monopole-antimonopole pair, and we 
have included a factor of (3^ to reflect both phase space and the velocity 
suppression of the magnetic coupling, as roughly implied by Eq. (5). 
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Fig. 3. Arrangement of the D0 tracking and transition radi 



Fig. 2: Arrangement of the DO tracking and transition radiation detectors. 

Any monopole produced at the Tevatron is trapped in the detector ele- 
ments with 100% probabihty due to interaction with the magnetic moments 
of the nuclei, based on the theory described in my review [9]. The experiment 
consists of running samples obtained from the old DO and CDF detectors 
through a superconducting induction detector. Figure 2 is a sketch of the 
DO detector. 

We use energy loss formula of Ahlen [15] to describe the interaction of 
the monopoles with the detector elements. 

Figure 3 is a diagram of the OU magnetic monopole induction detector. 
It is a cylindrical detector, with a warm bore of diameter 10 cm, surrounded 
by a cylindrical liquid N2 dewar, which insulated a liquid He dewar. The 
superconducting loop detectors were within the latter, concentric with the 
warm bore. Any current established in the loops was detected by a SQUID. 
The entire system was mechanically isolated from the building, and magnet- 
ically isolated by /i metal and superconducting lead shields. The magnetic 
field within the bore was reduced with the help of Helmholtz coils to about 
1% of the earth's field. Samples were pulled vertically through the warm 
bore with a computer-controlled stepper motor. Each traversal took about 
50 s; every sample run consisted of some 20 up and down t raver sals. Most 
samples were run more than once, and more than 660 samples of Be, Pb, and 
Al from both the old CDF and DO detectors were analyzed over a period of 
7 years. 
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Fig. 3: Sketch of the OU induction detector. Shown is a vertical cross section; 
it should be imagined as rotated about the vertical axis labelled "centerline." 

A monopole passing through the superconducting loop would produce a 
step in the current 



where L is the inductance of the loop, r is the radius of the loop, and a is 
the radius of the superconducting cylinder. The detector was calibrated with 
a pseudopole, a long solenoid, and the resulting steps in the output of the 
SQUID are seen in Fig. 4 to agree with theory. 

Figure 5 shows the histogram of steps from data collected from DO sam- 
ples. Similar histograms were obtained from the CDF data. 

From this, we can obtain limits on cross sections for the production 
of monopole-antimonopole pairs, and then, model dependent limits on 
monopole masses, as shown in Fig. 6. 

Table 1 shows the limits we obtained for different sample sets, and dif- 
ferent charges m', for various assumed production distributions. Our best 
mass limits are (assuming isotropic distribution) 




(15) 



• m' = \: fi> 265 GeV 

• m' = 1: n> 355 GeV 
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Fig. 4: Typical step plots: DO aluminum, CDF lead, and CDF aluminum. 
The experimental data was collected from pseudopole simulations; the steps 
shown are for the difference between the results with reversed polarizations 
of the pseudopole. Data agrees well with the SC theory, which incorporates 
the effect of the shielded superconducting loops. The theory without the 
shield, given by Barger and Ollson [16] is also shown. 



7 



60 




-4 -2 2 4 

Step (mV) 



Fig. 5: Steps from DO samples. A Dirac pole would appear as a step at 2.46 
mV. 



Table 1: Alternative interpretations for different production angular distri- 
butions of the monopoles, comparing 1 and 1 ± cos^ 6 (90% CL). Here the 
upper cross section limits corresponds to the distribution 1 + acos^ 9, and 



similarly for the lower 


mass limits n (all at 


90% confidence level). 
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Fig. 6: Cross section vs. mass limits. The three graphs show three different 
assumptions about the angular distribution, since even if we knew the spin 
of the monopole, we cannot at present predict the differential cross section. 
Shown in the second figure are the Bertani (1990) [13] and lunar (1973) [17] 
limits. 
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• = |: /!> 410 GeV 

• = 3: /i > 375 GeV. 



4 Reanalysis of Monopole Mass Limits Using Coulomb Re- 
summation 



We will now use the Solovtsov Coulomb threshold correction (4) in the form 

R = T{v)[Six)-\x{^ 
with T(v) given in (3) and S{x) given in (1), or 



R = T{v)[S{x)-lx{x)], (16) 



XU)=na^^^, Six) = Y^^. (17) 

We have simply neglected the perturbative term, as uncalculable. 
This is a small correction in QED, but here from Eq. (9) 

a^ml^^y, n= 1,2,3,.... (18) 

Figure 7 shows the substantial resulting increase in the cross section. This 
essentially pushes the cross section to the unitarity limit, 

^^.(2J + l)^3._ 

s s 

As a result, for all charge states, our limits become 

/x> 400 GeV. (20) 



5 Conclusions 

The relativistic Coulomb resummation factor plays an important role in anal- 
ysis of QCD experiments. Because the coupling is strong, it also plays a 
significant role in the theory of the production of magnetic monopole-anti- 
monopole pairs. Of course, because of the strong coupling, and even more 
because of the nonperturbative aspects of the Dirac string, there are poten- 
tially other effects which are just as strong but uncalculable. Our estimates 
of production rates were therefore extremely conservative, and a realistic as- 
sessment of the situation suggests that the limits on monopole masses from 
the Oklahoma experiment are at least as strong as the published limit from 
the very different CDF experiment [18]: 

/^CDF > 360 GeV, (21) 
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Fig. 7: Enhancement factor R for the theoretical monopole cross section due 
to Coulomb resummation, as a function of monopole mass = ^m, where 
A4 is the invariant mass of the monopole-antimonopole pair in TeV, and fi 
and m are in GeV. The bottom curve shows R given in Eq. (16), while in 
the top curve this factor is divided by the monopole velocity, since that is 
already included in the analysis, see Eq. (14). 
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